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Summary: The problem considered is that of the buckling of a right-angled isosceles 
triangular plate subjected to shear along the two perpendicular edges together with 
uniform compression in all directions. Four combinations of boundary conditions are 
considered, and the buckle is assumed to be symmetrical about the bisector of the right 
angle. In each case it is shown that, if the plate is subjected to shear only, the critical 
stress is changed considerably upon reversal of the shear. The larger value occurs 
when the shear is tending to decrease the right angle. Because of this fact, the inter- 
action curve, giving critical combinations of compression and shear, is unsymmetrical 
and the critical compressive stress can be appreciably increased by the application of a 

suitable amount of shear. 


Introduction 


In a recent paper"? the author considered the problem of the buckling of a 
clamped oblique plate (i.e. in the shape of a parallelogram) under pure shear 
parallel and perpendicular to two edges of the plate. It was found that the critical 
stress when the shear is tending to increase the obliquity of the plate is less than 
that corresponding to the reverse direction. It was also shown that for a 45° 
oblique plate the ratio between the two critical shear stresses lies between 1 and 
2-69 if the edges are clamped, and between | and 8 if the edges are simply-supported, 
the actual value depending on the side ratio of the parallelogram. 


A previous paper’) derived numerical results for the critical stresses of clamped 
oblique plates under compression parallel to two of the sides, and the question now 
arises as to the shape of the interaction curve, giving critical combinations of 
compression and shear. 


In a rectangular plate a reversal of the shear has no effect on its critical value, 
and the interaction curve is symmetrical about the axis of compressive stress. The 
curve is approximately parabolic, and the application of any shear stress reduces the 
critical compressive stress. For an oblique plate, however, the interaction curve is 
no longer symmetrical and it is therefore possible that the application of a certain 
amount of shear stress in a direction tending to decrease the obliquity of the plate 
may actually increase the critical compressive stress. 


The computation involved in the oblique plate problem is very tedious and in 
order to obtain some idea about the probable shape of the interaction curve a some- 
what easier problem, but having similar features, has been considered. This consists 
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45° 


Fig. 1. 


of a flat plate in the form of a right-angled isosceles triangle subjected to a 
combination of a shear stress g along the two perpendicular edges, together with a 
uniform compressive stress f (i.e. the same in all directions) as shown in Fig. 1. 
Four different sets of boundary conditions have been considered : — 


(i) All edges simply-supported, 

(ii) Perpendicular edges simply-supported, hypotenuse clamped, 
(iii) Perpendicular edges clamped, hypotenuse simply-supported, 
(iv) All edges clamped. 


In each case both the plate and the applied loading are symmetrical about the 
bisector of the right angle of the triangle. The buckling mode is therefore either 
symmetrical or anti-symmetrical about this bisector, and it has been assumed that 
the symmetrical type of buckle gives the lowest critical stresses in each case*. 


’ Sufficient numerical results were obtained to enable the interaction curve to be 
drawn for each of the first three sets of boundary conditions, as shown in Fig. 2. 
All three curves show similar effects; consequently no numerical calculations were 
performed on the plate with all edges clamped. 


It will be seen from Fig. 2 that the two critical shear stresses are considerably 
different in each case, and that the critical compressive stress is appreciably increased 
by the application of a suitable amount of positive shear stress (i.e. in a direction 
tending to decrease the right angle of the triangle). 


Although there is no mathematical analogy between the problem treated here 
and that of the oblique plate in combined compression and shear, nevertheless it is 
felt that it probably brings out the essential features of the interaction curve for an 
oblique plate. 


*This was proved by numerical calculation to be true for the case of all edges simply-supported. 
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Notation 
he 
ar x,y Cartesian co-ordinates (Fig. 1) 
lat a length of the equal sides of the triangle 
=x/a and y/a respectively 
ere Vv? = 07 /0€? + 07 
t thickness of plate 
oly E Young’s modulus 
ion 9 
D flexural rigidity = Er’ / {12 (1—+*)} - 
w lateral deflection of buckled plate 
t is f f uniform compressive stress applied to the plate 
— q__ Shear stress applied to the plate 
at , a’t 
ted. CS = =D f. =p4 respectively, buckling stress coefficients 
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m,n,p.,g integers 
=4(m—n), 4(m+n), 4(p—q), 4(p+gq) respectively 
Ann coefficient of w,,, in the series for w 


Wmn One of a,complete set of functions which satisfy the boundary 
conditions of the triangular plate 


Anpn, general element of the stability determinant 
Fe function defined by equations (16) 

Imp” definite integrals defined by equations (19) 
=litim=p, 
An, if (m+p) is even, =0 if (m+p) is odd 


2. Basic Theory 

Consider a thin plate, of uniform thickness ¢ and flexural rigidity D, in the 
shape of a right-angled isosceles triangle. Suppose that buckling has occurred 
under the action of a uniform compressive stress f in both the x and y directions and 
a uniform shear stress g, as shown in Fig. 1. The normal compressive stress acting 
on the hypotenuse is then equal to (f—gq), and there is no shear stress on this edge. 


If w is the lateral deflection at the point (x, y) and if every edge is either clamped 
or simply-supported, the condition that the potential energy should be a minimum 
is given by the equation 


{| ( DV ftV ?w+2qt bw dxdy=0 


where V ,?=0"/0x*+0?/dy’?, dw is any small variation in w which satisfies the 
boundary conditions, and the integration extends over the area of the plate. 


On transforming the variables x and y to the non-dimensional form € and 7 
by the equations 


f=x/a, 


where a is the length of the two equal sides of the triangle, the equation becomes 


¢ 
) Swdédy=0. . . 
0 0 
where V?=07/0€? + 0?/ On? 
a’t a’t 
and c= =D f. S= ‘ (2) 
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Suppose now that we represent the buckling mode by the series 


where the w,,, form a complete set of functions satisfying the boundary conditions 
of the plate. Then on substituting equation (3) into equation (1) and making 


we obtain an infinite set of homogeneous simultaneous equations in the a»; thus 


mn 


l 4 9 = 

V Wmnt CV? Winn S Woaq dédn. (6) 

A non-zero solution for the a@,,, exists only if the determinant formed by the 

coefficients Ampn, is zero. This gives a determinantal equation which must be 
satisfied by any critical combination of the buckling stress coefficients C and S. 


The problem now is to find suitable expressions for the functions w,, (€, ») 
‘| corresponding to each of the four sets of boundary conditions considered, and then 
to perform the integrations involved in equation (6) in order to obtain an expression 
for the general element of the stability determinant. Thereafter it is a matter of 
routine computation to obtain numerical results. 


3. All Edges Simply-supported 


To determine the functions w,,, in this case we notice first that the lateral 
deflection of a square plate, bounded by the lines €=0, 1 and »=0, 1, and simply- 
supported on all four edges can be represented quite generally by the series 


w (€, 1) = Sin sin ‘ (7) 
which satisfies the boundary conditions term by term. If this is to represent the 


deflection of the triangular plate it is clear that w must be anti-symmetrical about 
the hypotenuse =», 


i.e. w (f,n)+w £)=0. 
Hence (Ginn + Anum) Sin sin 
mn 
so that Qnm = — 
The series (7) now becomes 
) w (E, 1) (Sin sin N= — sin sin mr). 
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Numerical calculations showed that for all combinations of compression and 
shear considered, the mode giving the smallest critical stresses is symmetrical about 
the bisector of the right angle of the triangle. It is easy to show that this requirement 
is satisfied if the series contains only those terms in which (m+n) is odd. Also it 
is clear that interchanging m and n merely changes the sign of the general term of 
the series. Hence we have finally 


Sin sin Ney — sin sin ; (8) 

and we stipulate that m is even and n odd. 
On substituting equation (8), together with a similar expression for w,, (with p 
even and q odd), into equation (6), and performing the integration, we find that 


) if m=p and n=q 


2 2 2 
(m +n \(m +n C)+ — n? 2 
(9) 
(m? — q*)(n* — p*) 

Equation (8) is identical with that used by Klitchieff’, who considered the 
problem of the buckling of a right-angled isosceles triangular plate with simply- 
supported edges under pure shear stress. It will be shown later, however, that the 


numerical value which he obtained for the critical shear stress is not accurate. 


4. Perpendicular Edges Simply-supported, Hypotenuse Clamped 


In this case we again start with equation (7), giving a general expression for 
the lateral deflection of a square plate which is simply-supported on all four edges. 
In order that the series should represent the deflection of the triangular plate, it is 


evident that w must be symmetrical about the diagonal =» and also that w must 
be zero along this diagonal. 
From the condition of symmetry it follows that 
w (€,n)—w(n, &)=0. 
Hence, on using equation (7) 


> > Gini din) sin sin 0, 


so that 
Equation (7) now becomes 
w= Amn (sin sin + sin sin mrn) (10) 
mn 


where, in order to avoid repetition of terms, we stipulate that if (m+n) is even we 
make m=>n, and if (m+n) is odd we make m even and n odd. 
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It will be assumed in the numerical calculations that for all combinations of 
shear and uniform compression the governing mode is symmetrical about the 
bisector of the right angle of the triangle, so that 


w (&,n)=w(1—», 


It is a simple matter to show that this requirement is satisfied if every term in the 
series in equation (10) is such that (m+n) is even, and we therefore stipulate that 


We still have to ensure that w=0 along the clamped edge of the triangle, i.e. 
along =». Equation (10) then gives 


ann Sin sin 


mn 


so that Ann COS (m +n) TE= ay, COS (m—n) TE. 

On equating the coefficients of cosAré (where A=0,2,4....) in the series on the 
two sides of this equation, we obtain the following system of simultaneous 
equations : — 


It will be seen that the first of these equations can be disregarded, since it can 
be obtained merely by adding all the remaining equations together. We therefore 
solve the remaining equations to obtain expressions for the coefficients a,,,, in terms 
of the other coefficients of the type dmn(mAn). On substituting these expressions 
back into the series in equation (10) and collecting terms together we obtain finally 
the required expression for the deflection of the triangular plate; thus 


in which 

Wnn=Sin sin + sin nzE sin may + 2 sin 4 (m—n) sin 4 (m—n) 
(11) 

where (m+n) is even and m>n. 


It is a simple matter to check that the series satisfies term by term all the 
boundary conditions of the triangular plate. 
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On substituting equation (11), together with a similar expression for w,, (with 
(p+q) even and p>gq), into equation (6), and performing the integration, the 
general element of the stability determinant becomes 


32S mnpq mnpq 2 i?pq 
2r?mn 2s?mn 
2 j?s? ace arr ] 
+ Dir (i? + (7? sy i (7? ry (12) 
in which i=4(m—n), j=4(m+n), 
(13) 
r=4}(p-q), s=}(p+q), 
Terms of the type 4,,, and A,,, are defined by 
ari 
it 


Anp=1 if (m + p) is even, 


=0 if (m + p) is odd. . 


The numerical calculation of the elements of the stability determinant is not as 
arduous as it may appear from equation (12), since in every element a considerable 
number of terms are zero by virtue of the definitions of 4,,, and Any. 


Because of the restrictions on m, n, p and q it is clear that the parameters i, j, r 
and s are all positive integers. 


5. Perpendicular Edges Clamped, Hypotenuse Simply-supported 


In this case we notice first that the lateral deflection of a square plate, bounded 
by the lines €=0,1 and »=0,1, which is clamped on all four edges can be 
represented quite generally by the series 


w (6) F n(n) « 

where F,, (€)= sin sin ma 
(16) 

F, sin sin 


and every term satisfies the clamped boundary conditions. 
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If this is to represent the deflection of the triangular plate, it is clear that w 
must be anti-symmetrical about the diagonal €=y. It follows that dium—=—Qmn. 
and on assuming a mode which is symmetrical about the bisector of the right angle 
of the triangle, we obtain finally the expression 


War=F () F, —F, (é) Fe. () (17) 


in which m is even and n odd. 


On substituting equation (17), together with a similar expression for w,, (where 
p is even and q odd) into equation (6), the general element of the stability 
determinant becomes 


(18) 
1 
where Fe (6) F, () dé 
0 
1 


8 


0 


1 
Inj? =2 F,’ F, 
0 


and (d/d&) Fn (€), ete. 


The definite integrals defined by the equations (19) are given by 


=3 if m=p=1 
=2 if m=pA~1 
=-1 if m=p+2 
=0 otherwise. 
= —2(m’+1) if 


(m+p) if m=p+2 


otherwise. 


| 
= 
4 
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= 16 if me p= 1 
=2(m‘*+6m?*+1) if 


=— ag (m+ if m=p+2 


=0 otherwise. 


” 16 mp (m? + p? —2) 
(m? — p*) [(m? — — 8 (m? + p? — 2)] 


if (m+p) is odd 


=0 otherwise. 


6. All Edges Clamped 


In this final case we again start with equation (15), giving the deflection of a 
clamped square plate. The procedure is then exactly similar to that used in 
Section 4; thus we make w symmetrical about the diagonal =», and w=0 on this 
diagonal. We also assume that the buckling mode of the triangular plate is 
symmetrical about the bisector of the right angle. These requirements lead to the 
expression 


Wan =F» (€) (y) + Fa (©) Fn +2 F; (€) —2 F; (£)F; (y) : (20) 


where (m+n) is even and m>n, and i and j are defined by equation (13). 


On substituting equation (20), together with a similar expression for w,, (where 
(p+q) is even, and p>), into equation (6), the general element of the stability 
determinant becomes 


32S 


where 


+4 + Lie + 4 + Hie 
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— 2 Dig +2 +2 - 


In these expressions the definite integrals of the type /,,,)‘” are as defined by 
equations (19). Many of them will be zero, but in spite of this the task of setting 
up the stability determinant is still very tedious, and in fact has not been done 
beyond the second order. Once the determinant has been set up numerically, there 
is no more work involved in the calculation of the roots than for any of the other 
three sets of boundary conditions. 


7. Numerical Results 


For each of the first three sets of boundary conditions sufficient numerical 
results were obtained to enable an interaction curve to be drawn covering the full 
range of possible combinations of shear and uniform compression. 


The stability determinant is theoretically of infinite order but successive 
approximations to the roots can be obtained by including only a few appropriate 
rows and columns. As the order of the determinant increases the roots approach 
their exact values, the numerical magnitude being always too high. The roots were 
calculated by trial and error and the rows and columns which were included in the 
process of successive approximations were, in order, those corresponding to values 
of (m,n) and (p,q) given by the following sequences. 


(i) All edges simply-supported. 
(2, 1) (2, 3) (4, 1) (2,5) (4,3) (6, 1) (2,7) (4,5) (6, 3) (8, 1) (6, 5) 


(ii) Perpendicular edges simply-supported, hypotenuse clamped. _ 
(3,1) (4,2) (5, 1) (5, 3) (6,2) (7, 1) (6, 4) (7, 3) (8,2) (9, 1) ‘ 


(iii) Perpendicular edges clamped, hypotenuse simply-supported. | 
The same as (i). 


The numerical results are shown in Tables I, II and III, which also give the order 
of the determinant necessary to obtain sufficient accuracy, as indicated by the 
successive values of the roots as the order of the determinant increases. 


It will be seen that for the plate with all edges simply-supported the solution 
for pure compression (C=5, S=0) is obtained from a first order determinant. In 
fact this result is exact, for on substituting S=0 into equations (9) it is clear that 
every element of the stability determinant vanishes except those on the principal 
diagonal which have the value 


(m? + n*)(m? +n? —C). 


The value of the determinant is then equal to the product of the diagonal elements 
and one of these must vanish, giving C=m’+n’, where m and n must be chosen 
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TABLE I 
ALL EDGES SIMPLY-SUPPORTED 


| 
| 
| 


Compressive stress | 
coefficient C 0 | 5 


| 6:29 5-01 0 
Shear stress coefficient § — 11-55 0 11-57 18-51 30-89 
Order of determinant 8th Ist 9th 9th 11th 
TABLE II 


PERPENDICULAR EDGES SIMPLY-SUPPORTED, HYPOTENUSE CLAMPED 


Compressive stress 


coefficient C 0 7-82 8-23 5-02 0 
Shear stress coefficient § — 22-02 0 12-34 24-67 35-53 
Order of determinant 8th 4th 8th 10th 10th 

TABLE III 


PERPENDICULAR EDGES CLAMPED, HYPOTENUSE SIMPLY-SUPPORTED 


| 


Compressive stress 
coefficient C 0 9-35 10-98 7:55 0 


Shear stress = - 17-12 0 - 12-34 24-67 40-21 


8th Sth =| = 8th 10th =| 10th 


Order of determinant 


to give the smallest possible value of C. Since m is even and n odd, it follows that 
m=2 and n=1, so that C=5. The buckling mode then has the simple form 


w= sin sin sin sin2zy. 


A further result for the simply-supported plate, not shown in Table I, relates 
to the value of dC/dS at the point C=5, S=0. If we differentiate the stability 
determinant with respect to S and then put C=5, S=0, it follows immediately that 
at this point 


dC 128 
The results given in Tables I, II and III are shown plotted in Fig. 2, the broken 
line representing the tangent defined by equation (25). For all three sets of 
boundary conditions it will be seen that the two (positive and negative) critical 
shear stresses are considerably different, and that an appreciable increase in 
critical compressive stress is produced by the application of a suitable amount of 
positive shear. 

Because of the complexity of equations (21) to (24), the interaction curve for 
a plate with all edges clamped has not been derived. However, the point 
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corresponding to pure compression (i.e. S=0) has been calculated approximately, 
using a second order determinant, and is plotted in Fig. 2. The terms included in 
the series for w were those for which m=3, n=1 and m=4, n=2. Using the first 
term only gives C= 14-29, while both terms together give C=14:24. The closeness 
of these two values, and the experience with the other three sets of boundary 
conditions, suggest that this latter value is probably accurate to within one per cent. 


The problem of the buckling of a simply-supported triangular plate under pure 
shear has previously been considered by Klitchieff"’. He calculated only the 
positive critical shear stress and obtained the value S= + 46°32, as compared with 
S= + 30°89 given in Table I. The reason for this discrepancy lies in the fact that 
Klitchieff omitted several terms from the series for w which make a considerable 
contribution to the mode and he therefore obtained a value which is too high by 
about 50 per cent. 
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Subsonic Compressible Flow Past 
Bluff Bodies 


A. L. LONGHORN 


(Department of Mathematics, University of Manchester*) 


Summary: In this paper the Janzen“-Rayleigh”’ method is used to calculate the 
velocity potential for the steady subsonic flow of a compressible, inviscid fluid past a 
prolate spheroid. The fluid velocity at a point on the body is calculated. The analytic 
form obtained for this velocity differs, from that giving the velocity which an incom- 
pressible fluid would possess at the same point on the body, by a correction factor. The 
factor is an infinite series of first derivatives of Legendre functions of the first kind and 
odd order. The first three coefficients in this series are computed for bodies of certain 
axis ratios, and graphs of the values of these coefficients against axis ratio are plotted. 
The behaviour of the n‘" coefficient for large values of n is given. Results for slender 
ellipsoids, considering these as a limiting case of the family of ellipsoids just referred 
to, are obtained and are found to agree with the usual slender-body theory. Using 
these an attempt is made to continue the graphs of the first three coefficients in the 
correction factor series for the whole range of axis ratios of the ellipsoids in the 
system, namely zero to unity. The results obtained for the bluff-nosed ellipsoids may 
be used to estimate the effects of compressibility on the pressure distribution over the 

front of a general bluff-nosed body in steady flow. 


1. Introduction 


The Janzen"?-Rayleigh’’ method of calculating subsonic flows past obstacles 
by expanding the velocity potential in powers of the Mach number has been 
extensively applied to two-dimensional problems”). It has also been applied to 
three-dimensional flow by (a) Rayleigh®’, (6) Lamla™, (c) Schmieden® and (d) 
Sauer’. Rayleigh calculated the velocity potential for subsonic flow past a sphere 
and obtained a solution of the form ¢=9,+M’o,. Lamla has extended this work 
and has found a value for # including two further terms, those in M* and M’*. 
Schmieden has found the first two terms of the stream function for subsonic flow 
past an ellipsoid of revolution and has also calculated the velocity at the position 
of maximum thickness of the body and the critical Mach number of the flow. 
Sauer gives a numerical method which can be applied to calculate any axially 
symmetric three-dimensional flow. 


Our knowledge of the compressibility effects for flow past bodies of revolution 
is, therefore, limited to (i) the work already mentioned and (ii) slender bodies for 
which the usual perturbation theory can be applied. The present paper is an 
attempt to bridge the gap and in it the Janzen-Rayleigh method is applied to give 
some indication of the compressibility effects in uniform flow past bluff bodies of 
revolution. To formulate the problem mathematically a prolate spheroid is con- 
sidered as a typical bluff body and the flow past it is in a direction parallel to its 
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polar radius. By varying the eccentricity of the meridian section of this spheroid, 
nose shapes of varying curvatures are obtained. Thus if the nose region of any bluff 
body, for example an airship or an engine nacelle, can be approximated to in shape 
by a prolate spheroid, then the flow in that region of such a body can be estimated 
using the results given here. 


Only the first two terms in the expansion of the velocity potential in powers of 
the free stream Mach number M are calculated, the first term being the incompres- 
sible value of the velocity potential and the second the term in M*. However, for 
low values of M the neglect of M‘ and higher powers of M compared with M? gives 
a reasonable approximation. For flow at high Reynolds numbers the neglect of 
viscosity will lead to a solution which gives a reasonable picture of the flow over the 
front portion of the body, for here the actual flow approximates fairly closely to 
that calculated from potential theory. In practice the effects of viscosity are 
manifest mainly towards the rear of the body where the boundary layer begins 
to thicken and finally separate, thus producing the wake. Hence, on the rear of 
the body the solution obtained by neglecting viscosity cannot be expected to give 
any useful predictions as in this region the viscous effects are the dominant feature 
of the flow. However, the separation at the rear is known not to affect greatly the 
pressure distribution over the front of the body and in view of this fact the solution 
found in this paper may be expected to give a reasonably accurate representation 
of the compressibility effects in the front of the body. 


In the application of the Janzen-Rayleigh method to this particular problem 
the equation satisfied by the velocity potential # in steady subsonic flow, namely 


2 
(where a is the speed of sound through the fluid at a point where the density of the 
fluid is p and its speed is q) is expressed in ellipsoidal co-ordinates ((,“,). The 
fluid is assumed to obey the adiabatic law, and hence 


where U is the free stream velocity, a,, is the speed of sound through fluid moving 
with this velocity and of density p., and y is the ratio of the specific heat of the 
fluid at constant pressure to that at constant volume. For the flow considered, 
the free stream velocity is U in the negative Ox-direction and parallel to the polar 
axis of the ellipsoid and hence M, the free stream Mach number, is given by 
M=U/a,,.. If now, a series expansion for ¢ in terms of M is assumed, then 


where 9, is the value of ¢ when M=0, that is the solution for the incompressible 
flow of a stream with velocity U past the given body. In this paper, only 9, is 
calculated so that the value of '» which is found has the form 
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The expression obtained for 9, is of the form 


where the coefficients b.,,, are independent of », and P24, (“) is the Legendre 
coefficient of the first kind. The fluid velocity at a point on the body is then 
calculated and it is expressed in terms of the velocity which the fluid would possess 
if it were supposed incompressible. These results, together with some computations, 
are described and discussed fully in Section 6. 


Notation 
A 


A.R. 
(C,)o = 
(Cy)o. «1 a 


(C,)m. Bya 


F - 


F, (1) 


F,Q.F:Q.FsO 


Ans; 
M = 

Pans 
Pons’ = 

(Y = 

U 

a 

ax 


constant; see equation (8) 
axis ratio 
(p— px.)/(4 pressure coefficient for incompressible flow 


pressure coefficient at a point on a body of thickness ratio 
B(1—M?)!/ in incompressible flow 


pressure coefficient at a point on a body of thickness ratio 
8/z in compressible flow 


t?/2 in case of slender bodies; see equations (49) and (50) 
[Hogi (Gy) Pong’ case of bluff ellipsoids; see 


equation (29) 
F («) in case of slender bodies 
Gy) in case of bluff ellipsoids 
n=0 
see equation (11) 
see equation (12) 
see equations (20) and (22) 
see equation (28) 
U/a,,, free stream Mach number 


Legendre coefficient of first kind 

d 

du ones 

Legendre coefficient of second kind 


d 


free stream velocity 
speed of sound at any point in fluid 
speed of sound in free stream 
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polar radius of spheroid 

see equation (5) 

equatorial radius of spheroid 
constant; see equation (24) 
constant 


ratio of specific heats of fluid at constant pressure and 
constant volume 


constants 

see Section 5 

eccentricity of meridian section of spheroid 

see equations (12), (13), (14) and (16) 

velocity potential for compressible flow past the spheroid 
velocity potential for incompressible flow past the spheroid 
correction terms to be added to ¢, to obtain @ 

perturbation potential; see Section 5 

see equation (26) 

constant 


angle between the tangent to the meridian section of the 
ellipsoid and the x-axis 


static pressure at any point in the fluid 

Static pressure in free stream 

speed of fluid at any point in compressible flow 
velocity of fluid at any point in compressible flow 
speed of fluid at any point in incompressible flow 
velocity of fluid at any point in incompressible flow 


speed of fluid at point on body surface in incompressible 
flow 


speed of fluid at point on body surface in compressible flow 
radius of section of spheroid perpendicular to x-axis 
density of fluid at any point 

density of fluid in free stream 


length of arc of a meridian section of the body measured from 
the forward stagnation point; see equation (52) 


thickness ratio 
components of perturbation velocity; see Section 5 
Cartesian co-ordinates 


ellipsoidal co-ordinates; see equation (6) for relation between 
these and (x, y, z) 


a particular value of ¢, defining a prolate spheroid 
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2. Derivation of the Velocity Potential of the Flow 


Consider a prolate spheroid with polar and equatorial radii « and £ respectively. 
Choose rectangular Cartesian axes Oxyz at O, the centre of the spheroid, and fixed 
in space, and let the polar radius coincide with Ox. Suppose a steady stream with 
velocity U in the negative Ox-direction flows past the body. With co-ordinates 


x=Kp, z=K 1) sinw, (6) 


where K is a constant and (x, y, z) are the Cartesian co-ordinates of a point referred 
to the axes Oxyz, then this point may also be denoted by (¢,,«), the so-called 
ellipsoidal co-ordinates. These latter form an orthogonal system and in particular 
(={,, a constant, denotes an ellipsoid of revolution with «=K(, and B=K ((,?— 1). 
In terms of the co-ordinates ((,,) the expression obtained for 9%,, the velocity 
potential for the motion of a steady stream of incompressible fluid with velocity 
U in the negative Ox-direction past the ellipsoid of revolution (=(,, is given by” 


where the fluid velocity q,—grad¢, (the flow being irrotational), Q,(¢) is the 
Legendre coefficient of the second kind and A is a constant, namely 


where e is the eccentricity of the meridian section of the ellipsoid (=1/(,). 


For a compressible fluid in steady motion the flow is irrotational and hence a 
velocity potential exists. satisfies the equation (1) 


1 
grad grad q’, 
and since in ellipsoidal co-ordinates the operator V? is given by 


and here the motion is symmetrical about Ox, and thus ¢ is independent of », then 
equation (1) becomes 


The boundary conditions for the flow are 


(i) the normal fluid velocity on the surface of the ellipsoid is zero, that is 
09/0¢=0 when (= , for all values of », and 


(ii) at large distances » tends to the value of the velocity potential for a free 


stream flowing with velocity U parallel to the negative Ox-direction, that is 
¢—> —UK for large values of ¢. 
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Assume a solution of equation (9) of the form ¢=9¢,+M’9,, where M(=U/a,,) is 
the free stream Mach number. From equation (2) 


2a? + 1) (y— 1)U?, 


1 M? 


|, 


and hence, if in equation (9) the expressions for 1/a* and @ are substituted, then 
equating coefficients of M? gives 


since V*¢,=0, where q, is calculated using the velocity potential ¢,. Evaluating the 
right hand side of equation (10) gives 


(11) 
where F,(0)=3{UK6+AQ, 42, O}. 
and F, (= {UK(+ AQ, (0)} fe 1) {UK(+ AQ, 17 


where f, (Q and f ,(Q) are functions to be determined and substitute in equation (11). 
If f, (©) is chosen to be a solution of the equation 


— _ Ay} 
and f, (0) is chosen as 
2 A? 
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then the equation for 2, becomes 


+5 } = let -20, (0. 


(15) 
From equation (13), 
+40,’ OF dé-O,@ | P, OUK + AQ, OF dé. 
co 
(16) 
To solve equation (15) consider the expansion 


which follows from the expansion given in Ref. 8. For this problem the expansion 
is certainly valid, as here ¢ is real and greater than unity for all surfaces of the 
family except the degenerate case, which is of no interest, and also » lies between 
—land +1. Differentiating (17) with respect to ¢, 

3) 


Qonai Pansi . (18) 
Hence ‘igi in equation (15), 


0 { 


AP 4n+3 20, (0) 


Assume a solution for ?, of the form 


4n43 


where G.,,,, (() are functions to be determined, and substitute this series in equation 
(19). Equating coefficients of P2,4;(«), we then obtain 


~{) — (2n+ 1) 2) Gang (Q= ] 


. ~ 
Hence 
(0) Oona (2) {a 2n+1 a¢3) } dé — 

(22) 


this being true for n=0, 1, 2,... 
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As a particular integral for equation (11), we then have 


A 
where the functions f, (9), f. (G) and G.n4,(@ (n=0, 1, 2, . . .) are given by equations 
(16), (14) and (22) respectively. The general solution of equation (10) can then be 
written as 


the last two terms being the complementary function for equation (9) and c,, d,’ 
being constants to be determined by applying the boundary conditions (i) and (ii). 
Now ¢,, being the solution for incompressible flow, will satisfy both (i) and (ii) and 
thus ¢, must also separately satisfy (i) and must tend to zero at large distances from 
the ellipsoid, that is as ¢ tends to infinity. Hence, from the latter condition 


n=0 


must tend to zero as ( tends to infinity, since all the terms of the particular integral 
do so. This can only be true if c,=0 for all nm. The constants d,’ still remain to be 
determined, and for this purpose (i) must be used. We have 


A? n 
+ = (4n + 3) Pons (| } t=t + 
Equating coefficients of P,,(u), then d,’=0 for all even values of n. Writing 


A 4n+3 d 
a,2K? 2 2n+1 


Sa? ad 


= 


then we have, for n=0,1,2...., 


+3 +3 d 
Hence 
and (25) 
d ] A? (Q 1 ] 
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for n=1, 2,3... . These relations determine the constants d,,,, and hence ¢ is 
determined. The form of the solution is 


A A 


x { fo Qon4i (Q + + 0 Q } 


n=0 
where C.n+,’ are constants and g,n+, (¢) are known functions of ¢. 


3. Calculation of the Fluid Velocity at a Point on the Body 


To obtain the fluid velocity at a point on the body only (0¢/0,);-;, is required, 
since (69/0) ;-;, is zero from the boundary condition (i). From equation (26) 


og = 20 A A AP on 41 (1) 
A* don 
Gantt (<)+ Qon+1 © | 
Hence 
0%, 


co 


A 
» 


having substituted for G,,+, (¢,) from equation ed Simplifying the last expression 
still further by substituting for f, (¢,) and f, (¢,), then 


= (4n+ (Go) (60) + |, 


= (2) A 
Op t= Op. a,,*K? 


% 
1)0,’(,) | QO, (©) [UK + AQ,’ 


co 


% 
J Ones (€) Ome ae } 


But = —A /(&2—1) and hence the expression for (09/0);-;, becomes 
after simplying 


, 
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Fig. 1. 


The functions H, ({,), H, (¢,) and H, (¢,) plotted against axis ratio. 


% 
—3 0,7 (6) + 1) 2 30,° 


, 12 4 3 5 3 0 42 
+Q,7Q, 3 > 72: + +Q,Q, ] 


(28) 
and 


% 
(4n 3) 1 Qon41” 1) = Qon41” 
Any, G)= 4 { | —1) dé D 


forn=1,2... 


The value of qv, », the fluid velocity at a point on the surface of the body, is 


_ 
K — py Op 


Hence if q,,, is the velocity of the fluid at the same point on the body when the 
fluid is supposed incompressible, then from equation (27) we have 


a= 


qv, m= Yb, 0 {1 +M? 2. { Hons: ()} ]. (29) 
Values of H, ((,), H; ((,) and H, (¢,) for varying ¢, have been computed and these are 
given in Table I and also plotted in Fig. 1. 
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TABLE I 
the body So H, H; (G) 
1 ore) 0-117 —0-061 0 
0-81 1-70 0-092 — 0-044 —0-0039 
0-60 1:40 0-065 0-027 — 0-0056 
0-42 1-10 0-040 —0-014 — 0:0048 
0-20 1-05 0-013 —0-0027 —0-0016 


4. Behaviour of H+ (%) for Large Values of n 


It is known? that exp (— x) 


for large values of n, where x=n(u?—1)! and »> 1. Hence, applying this result we 
have for large n 


(2n ra «,? 1)! exp [ (2n 1) 1)*] (30) 


where {, is real, fixed and greater than unity. 


oo 
(E)(& + I) dé 

Consider 


and denote this by I,,4,(¢,). 


en anti \So 2(2n+1) 


exp [—2(2n+ 1) 1)4] d€=1(,) 


2 (2n+1) - (2-1) exp[—2 (2n+ 1) — dé 
1) 


~ 4(Qn4+ 1? exp[—2(2n+ 1) — 1)#]. 


+1) 


Hence | ~I 4 (2n rs 1)? a, = 


exp [- 2(2n+1)((,?- 14]. GBD 


Further, we have 


exp[—2(2n+1)@?-D4] G2) 


and hence from (31) and (32) it is seen that for large values of m the term in the 
expression for H,, (¢,) which involves the integral may be neglected compared with 
the other term in the expression, these terms being at least in the ratio of 1/n: 1. 
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Differentiating (30) with respect to ¢,, we have, for large values of 7, 


(2n+ 


Hence, using the results of (31), (32) and (33), it follows that 


The coefficients H,,,,, (¢,) are thus seen to behave like n-! exp[—(2n+ 
for large values of n. 


5. Thin Ellipsoids of Revolution 
For the incompressible flow of a steady stream of fluid with velocity U past the 
ellipsoid of revolution (={,, the velocity potential », is given by” 
UKn2, 


The perturbation potential ¢, on the steady stream flow due to the presence of the 
ellipsoid is therefore given by 


Py 


_ UKp2, 


If (u,, v,) are the components of the perturbation velocity in the Ox and radial 
directions respectively then 


_ 


osa, 


where A is the angle that the tangent to a meridian section of the body makes with 
the Ox axis. Clearly by using (6) it follows that cos A= (1 — and 
hence 


For a slender body ¢, is approximately unity (for example for a body of axis ratio 
1: 10, §,=1-005). Hence write ¢,=1+8 where 8 is small. Substituting in (37), 


1+5, 2+8 


(1 — p? + 26 + 8?) 
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For slender bodies the axis ratio 8/« is of the order of 0-1 or less and hence 8 is less 
than or equal to 0-005. Thus 


8 < 
1-p? 


if 1—p»?>0°05, that is if »«=<0-97. 


Hence certainly for » <0-95 we may expand in powers of 6/(1—,*) and the result 
will be sufficiently accurate when we neglect squares and higher powers of this 
quantity. Thus, neglecting 5? compared with 6 in the denominator of (38), and then 
expanding the denominator in powers of 25/(1—?), we have 


Since (,=1/e=1+4, then Thus and, 
on writing t=8/2, equation (39) becomes 


2 
If the equation of a section of the ellipsoid is x*/a?+r,?/8?=1, then since x=ap, 


we have But 1—y»?>0-05 and therefore, r,?/8?=>0-05. Hence 
t?/(1—”) is at most of order t. Equation (40) then becomes 


2 3 


The boundary condition on the ellipsoid is 


=-tana= (1-2). 


But u,=O (U?’ log t) and hence 


since (t8/r,)(1—r,?/8)' is at most of order r'. If g,,, is the speed of the fluid at a 
point on the body surface, then from (41) and (43) 


b o =U 


=u" [ {1-1 (1-10g2) +0 {2 (1- } 


=U? [ (1 =) + (1- log) |. 


2 
Hence 1- = log +e (3- log . (44) 
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If (C,), is the pressure coefficient and is defined as 


Che 
Pa" 


where p,, is the pressure in the free stream, then 


Using the result for slender bodies") that 


(1—-M?*) 
which relates the compressible flow past a body of axis ratio 8/a with the incom- 
pressible flow past a body of axis ratio B(1 —M*)#/2, and (C,) 
being the pressure coefficients referred to these two flows, then we have from 
equation (44) the result 


2 2(3_ & 


where q»,™ is the speed of the fluid at a point on the slender body of axis ratio B/a 
when the fluid is compressible. From (44) and (47) it follows at once that 


+28? log {2/[t (1 — — /r,”) 
1 + 22? log (3 — B?/r,”) 
=1-f 


Hence we have 
Qv,o” [1 — log (1 — M’)], 


and expanding in powers of M? and neglecting M* and higher powers we have, to 
the required degree of approximation, 


Comparing this result with the general form of the result for bluff bodies, namely, 


it is seen that for slender ellipsoids F() is a constant, in fact F(u)=f?/2= 
B?/(2a?)=8. It must be remembered that this result has been shown only for the 
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range » <0°95, that is for x < 0-952. For points nearer the end of the body the 
slope of the body section increases rapidly and finally, at x=a, becomes infinite. 
This factor renders invalid the method here adopted, as it is necessary for the slope 
of the body surface to be small for the method to hold. The result obtained in 
equation (49) is easily obtained by the usual slender-body theory using source-sink 
distributions on the body axis. 


6. Discussion of Results 


It has already been stated in the Introduction that the form of the solution found 
here is 


¢=%,+M’o, where = anes 


ju=x/2, the coefficients b,,,, which are independent of » and P»,,,() being the 
Legendre coefficients of the first kind. The fluid velocity g,,™ at a point S on the 
surface of the body is calculated and it is given by 


dv, M =dp, 0 [1 M? { (Go) dp Pandy (u) } ] (51) 
where (¢, is the parameter denoting the particular spheroid considered (this being one 
of a family of confocal prolate spheroids) and q»,, is the velocity which the fluid at 
S would possess assuming it to be incompressible. Hence 


d 

n=0 d 
is a correction factor by which q,,, must be multiplied in order to take into account 
the effect of compressibility. The values of the coefficients H, (¢,), H, (¢,) and H, (4) 
occurring in equation (50) have been calculated for some values of ¢, (that is for 
bodies of various axis ratios) and are plotted against axis ratio. These results are 
given later. The behaviour of H.,,,, (¢,) for large n is given, and this coefficient is 
found to decrease exponentially. 


The parameter (, in the co-ordinates here used defines a prolate spheroid. For 
the whole family of confocal prolate spheroids ¢, varies in value between unity and 
infinity. When ¢,=1 the ellipsoid degenerates into a straight line joining the foci of 
the family, and when (, tends to infinity the ellipsoid tends to a sphere of infinite 
radius. However, the case of (,=1 is clearly of no interest and thus in all cases 
considered here ¢, will be greater than one. Furthermore, the case of a sphere has 
been discussed by Rayleigh’’) and hence the results when (, tends to infinity can be 
obtained without using the present method. For an ellipsoid of axis ratio 4: 5, (, 
is only of the order of 1:7, so no large values of (, need enter into the computations. 
As (, approaches unity, computation becomes more laborious as the integrals 
required have a singularity at (,=1. Even for bodies of axis ratio 1 : 5, when 
{,=1-02, the calculations are quite tedious and this factor limits the number of 
results which are reasonably obtainable for bodies of smaller axis ratio. Further- 
more, as the axis ratio decreases it becomes necessary to calculate more coefficients 
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BLUFF BODIES 


TABLE Il 
VALUES OF THE FUNCTIONS PLOTTED IN FIGS. 2 AND 3 


u=1 or s/a=0 give the stagnation point on the nose of the body, the arc length s 
being measured from this point 


Sphere 1 0-9 0-8 0-6 0:4 0:2 0 
Axis Ratio| s/2 0 0:4538} 0-6458 | 0-9258 | 1-1606 | 13701 | 1-5708 
= 1-00 | —0:249 | —0-1621 | —0-0843 | 0-0438 | 0-1353 | 0-1902 | 0-2085 
Ellipsoid 0-9 0:8 0-6 0-4 0-2 0 
Axis Ratio| s/2 0 0:3732| 0-5393 | 0-7931 | 10182 | 1-2239| 1-4245 
=0°81 F,(v) | —0-2305 | —0-1344 | —0-0579 | 0-0488 | 0-1103 | 0-1413 | 0-1507 
Axis Ratio | 0-9 0:8 0:6 0-4 0:2 0 
=0-60 0 0:2854| 0-4287 | 0-6590 | 08755 | 1-0770 | 1-2774 
F,(v) | —0-481 | —0-0946 | —0-0308 | 0-:0464 | 0-0805 | 0-0924 | 0-0950 
AxisRatio| yp 1 0-9 0:8 0-4 0-2 0 
=0-42 s/a 0 0:2149| 0:3360| 0:5495 | 0-7586 | 0-9591 | 1-1593 
F,(v) | —0-116 | —0-0560 | —0-0120 | 0-0351 | 0-0505 | 0-0525 | 0-0520 
Axis Ratio| 1 0-9 0:8 0-6 0-4 0:2 0 
=0:20 0 0:1439| 0:2463 | 0-4473 | 0-6504 | 0-8510 | 1-0510 
F,(») | —9-0271 | —0:0095;| 0-0022 | 0-0138 | 0-0160 | 0-0149 | 0-0141 
Axis Ratio m ] — — — — — 0 
F, (2) is a constant and has the value 0-005 | 


in the series given in (50), for the further coefficients become increasingly comparable 
with the first three. Thus, for bodies of axis ratio less than about 2:5 more 
coefficients than three should be found and, in addition, the calculation of each is 
worsening, so that the results given will be useful chiefly for bodies of axis ratio 
greater than 0-4. Accordingly, the functions H, (¢,), H,(¢,) and A, (¢,) as given by 
equation (28) have been calculated for values of ¢,=1-05, 1:10, 1-40 and 1-70. 
These values correspond to bodies of axis ratio 0-20, 0:42, 0-60 and 0-81 respectively, 
all these figures being given correct to two decimal places. These results, together 
with the case of the sphere'’’ are given in Table I and plotted in Fig. 1. 


The case of the ellipsoid becoming slender is considered as the limiting form 
of the solution obtained in (50). (, is replaced by 1+4 where 4 is small and so 
H, H; and are obtained for bodies of axis ratio 0-1 and less. Com- 
bining these results with the ones already found and using Rayleigh’s result for 
a sphere’, complete curves of these three coefficients for bodies of axis ratio from 
0 to 1 can be drawn. A reasonably smooth join is obtained between the portion of 
the graph representing the slender-body results and that giving the bluff-body 
results (Fig. 1). 


If the general form of the velocity q,,» at a point S on the body is 
=v, o [1+ (»)], 
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Fig. 2. 
The function F,, («) plotted against » for bodies of different axis ratio (A.R.). »=1 gives 
the stagnation point at the nose of the body. 


then Fig. 2 shows the function F,(u) plotted against » for bodies of certain axis 
ratios, where F,, («) is an approximation to F (u), namely, the first three terms in the 
series for F(u). The results are also displayed in Fig. 3 by plotting F, (u) against 
s/a, where s is the arc length of a meridian section of the body measured from the 
forward stagnation point and 2 is the semi-major axis of the body; s is, of course, 
given by an elliptic integral. We have 


and hence s= =| ao, 
1 
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-0:25 ‘ 
Fig. 3. 
The function F, («) plotted against s/a for bodies of different axis ratio (A.R.). 


Using Fig. 1 it is possible by interpolation to read off the values of the 
coefficients H, (¢,), H, (¢,) and H, (¢,) for a given body, that is, given the axis ratio 
of the ellipsoid which best approximates to the body, and then to plot F, () for that 
body. Alternatively, it is possible to interpolate for F, (u) from the graphs of this 
function which are given. As an example of the latter method, consider the case of 
a body of axis ratio 0-5. The value of F,(u) will be approximately —0-15 when 
s/z=0, it will be 0 when s/2=0-45 and will rise to +0-07 when s/a attains its 
maximum value for this body. 
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